Bulk effects in the cosmological dynamics of brane-world scenarios 
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In this paper we deal with the cosmological dynamics of Randall- Sundrum brane-world type scenar- 
ios in which the five-dimensional Weyl tensor has a non-vanishing projection onto the three-brane 
where matter fields are confined. Using dynamical systems techniques, we study how the state 
space of Friedmann-Lemaitre-Robertson- Walker (FLRW) and Bianchi type I cosmological models 
is affected by the bulk Weyl tensor, focusing in the differences that appear with respect to general 
relativity and also Randall- Sundrum cosmological scenarios without the Weyl tensor contribution. 
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I. INTRODUCTION 

Recently, Randall and Sundrum have shown that for 
non-factorizable geometries in five dimensions the zero- 
mode of the Kaluza-Klein dimensional reduction can 
be localized in a four-dimensional submanifold where 
Newtonian gravity is effectively reproduced at large dis- 
tances [jjj. The picture coming out is one with all mat- 
ter and gauge fields, excepting gravity, confined in a 
three-brane embedded in a five-dimensional spacetime 
with Z2 symmetry. This scenario was motivated by orb- 
ifold compactification of higher dimensional string theo- 
ries. In particular by the dimensional reduction of eleven- 
dimensional supergravity in R 10 x S 1 /^ introduced by 
Hofava and Witten ||. The cosmological implications 
of the Hofava- Witten theory have been analyzed in || . 
The actual relevance of this type of scenarios is that they 
provide an alternative framework to investigate new so- 
lutions to the hierarchy and cosmological constant prob- 
lems. 

Gravity on the brane can be described by the Ein- 
stein equations modified by two additional terms, one 
quadratic in the matter variables and the second one be- 
ing the "electric" part of the five-dimensional Weyl ten- 
sor 0-0. In a previous work || we studied the whole 
dynamics of the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) and the Bianchi type I and V cosmological mod- 
els taking into account only the effects produced by the 
first type of corrections. In that paper we found new 
critical points corresponding to the Binetruy-Deffayet- 
Langlois (BDL) models Q, representing the dynamics at 
very high energies, where the extra-dimension effects be- 
come dominant. The presence of these models seems to 
be a generic feature of the state space of more general cos- 
mological models. Moreover, we saw that the state space 
presents new bifurcations and how the dynamical charac- 
ter of some of the critical points changes with respect to 
the general-relativistic case. Finally, we showed that for 
models satisfying all the ordinary energy conditions and 
causality requirements the anisotropy is negligible near 



the initial singularity. On the other hand, for late times 
the anisotropy is diluted as in general relativity (GR) 
but now there are intermediate stages in which the shear 
contribution grows. 

Here, we will complete the analysis done in || (here- 
after Paper I) studying the effects of the five-dimensional 
Weyl tensor on FLRW and Bianchi type I cosmological 
models. In the case of FLRW models this study will be 
completely general whereas in the Bianchi type I case we 
will neglect the Weyl tensor components for which the 
theory does not provide evolution equations. The pa- 
per is organized as follows. In Sec. || we briefly describe 
the geometric formulation of brane-world scenarios and 
the dynamical equations for FLRW and Bianchi type I 
cosmological models with a perfect-fluid type matter con- 
tent. The dynamics of the FLRW and Bianchi type I cos- 
mological models will be analyzed in Sees. Ill and IV re- 



spectively. We will finish with some concluding remarks 
111 Sec. 0. Through this paper we will follow the follow- 
ing notation: upper-case Latin letters denote coordinate 
indices in the bulk spacetime (A, £?,... = — 4) whereas 
lower-case Latin letters denote coordinate indices in the 
brane (a, 6, . . . = — 3). We use physical units in which 
c = 1. 



II. PRELIMINARIES 

In what follows we introduce the geometric framework 
to analyze the brane-world scenario and the main as- 
sumptions we will adopt to study the FLRW and Bianchi 
type I cosmological models. 



A. Geometric formulation of brane-world scenarios 

In Randall-Sundrum brane-world type scenarios mat- 
ter fields are confined in a three-brane embedded in a 
five-dimensional spacetime (bulk). It is assumed that 
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the metric of this spacetime, g\ Bl obeys the Einstein 
equations with a negative cosmological constant A( 5 ) 
(see §01) 
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where G AB denotes the Einstein tensor and K(k\ is the 
five-dimensional gravitational coupling constant. More- 
over, Tab is the matter energy- momentum tensor and 
the Dirac delta function reflects the fact that matter is 
confined in the spacelike hypersurface x A = x — 0, the 
three-brane, with induced metric qab an d tension A . 

From the Gauss-Codacci equations, the Israel junc- 
tion conditions and the Z2 symmetry with respect to the 
three-brane, we can see §0] that the effective Einstein 
equations on the brane are 



Gab = -Agafc + K 2 T ab + K( 5 )£?o6 
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where G ao is the Einstein tensor of the induced metric 
g ao . The four-dimensional gravitational constant n and 
the cosmological constant A are given in terms of the 
fundamental constants in the bulk (k( 5 ),A( 5 )) and the 
brane tension A. It is important to note that in order 
to recover conventional gravity on the brane A must be 
assumed to be positive. As one can observe, there are 
two corrections to the general-relativistic equations. The 
first one, denoted by the term S a b, is quadratic in the 
matter energy-momentum tensor T ao (T = g ab T ao ) 
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±g ab [3T cd T cd - T 2 ] . (3) 



The second correction, E ab , correspond to the "electric" 

part of the five-dimensional Weyl tensor, C ABCD , with 
respect to the normal, ha (n A nA = 1), to the hypersur- 
face x = 0, that is 



p (5) _ r (5 
^AB — ^ACBD 



C D 

n n 



From the modified Einstein equations (g) and the energy- 
momentum tensor conservation equations (V a T a b = 0) 

we get a constraint on S a b and E a 
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V a (4?-«(6)^) = 0. 
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Following |9|], it is useful to decompose E^ b with re- 
spect to any timelike observers u (u a u a = — I) in a scalar 

part, U, a vector part, Q a , and a tensorial part V a b, in 
(5) 

such a way that E ab can be written as 



E 



(5) 



(5) 
K 



[(u a U b + \h ab )U + 2U( a Qb) + Vab] 



(5) 



where the following properties hold 

Q a U a = , V { ab) - Vab , V a a = , V ab U b = . 



That is, Q a is a spatial vector and Vab is a spatial, sym- 
metric and trace-free tensor. The scalar term has the 
same form as the energy-momentum tensor of a radia- 
tion perfect fluid, and for this reason U was called "dark" 
energy density in IJ. Note that the constraint (|J) pro- 
vides evolution equations for U and Q a , but not for V a b 
(see [§). 



B. Cosmological dynamics on the brane 

Throughout this paper we will describe matter by a 
perfect fluid 

Tab = (p + p)u a u b + pg ab , 

where u, p and p are the unit fluid velocity of matter, the 
energy density and the pressure of the matter fluid re- 
spectively. We will also assume a linear barotropic equa- 
tion of state for the fluid 



P= (7- !)/£>• 



(6) 



The weak energy condition (see, e.g., |lfj]) imposes the 
restriction p > and from causality requirements we have 
that the parameter 76 [0, 2]. 

In this paper we deal with non-tilted homogeneous cos- 
mological models on the brane, i.e. we are assuming that 
the fluid velocity is orthogonal to the hypersurfaces of 
homogeneity. Then, it is convenient to make the split- 
ting (g|) with respect to the fluid velocity. In particular 
we will consider: (i) The FLRW models, with metric ten- 
sor given by 



ds^ 



-dt 2 + a 2 (t) [dr 2 + E 2 k (r)(d9 2 + sir 



0] 



where a(i) is the scale factor, k is the curvature param- 
eter, and 

sin r for k = 1 , 
E fc (r) = ^ r for k = 0, 
sinh r for k = — 1 ; 

and (ii) Bianchi type I models, whose line element can 
be written as 



ds^ 



^dt 1 
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Then, in both cases the fluid velocity is u = d/dt. We 
will not consider Bianchi type V models as we did in Pa- 
per I because, as we showed there, their main dynamical 
features are a combination of those of the FLRW and 
Bianchi type I models. 

Taking into account these assumptions and the effec- 
tive Einstein's equations (^) , the consequences of having 
a FLRW model on the brane are 

Qa = Vab = , 
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and this, through the constraint (Q), further implies 
D a U = <S> U=U(t), 



(8) 



where D a denotes the covariant derivative associated 
with the induced metric on the hypersurfaces of homo- 
geneity [h a b = g a b+u a Ub). Instead, in the case of Bianchi 
type I we get Q a — but we do not get any restriction 
on Vab ■ Since there is no way of fixing the dynamics of 
this tensor we will study the particular case in which it 
is zero. Then, also the condition (||) follows. 

The key equation to describe the dynamics of these 
models is the Friedmann equation, which determines the 
expansion of the universe, or in other words, the Hubble 
function: ZH = V a u a . The general expression of this 
equation in the case of homogeneous cosmological models 
is 
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where 3 R is the scalar curvature of the hypersurfaces or- 
thogonal to the fluid flow and 2a 2 = a ab a a b is the shear 
scalar (<r a b = h a c hb d S7 ( c u d j — Hh a b). We consider only 
the case of a positive fourd-dimcnsional cosmological con- 
stant, i.e. A > 0. For Bianchi type I models 3 R vanishes 
whereas for FLRW models it is given by 3 R = 6ka~ 2 (t) . 
On the other hand, the shear vanishes for FLRW models 
and for Bianchi type I models the evolution of the shear 
scalar is 

(a 2 )' = -6Ha 2 . 
In both cases the evolution equation for U is || 
U = -4HU , 

and, from the energy-momentum conservation equation, 
the evolution of p is simply 
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Finally, an important quantity in the dynamical study 
will be the deceleration parameter q, defined by 



H = -{l+q)H 2 



(11) 



where H is found to be 



H = -H 2 -^n 2 p 
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which is the modified Raychaudhuri equation. 



effects of the extra dimensions are only due to the term 
quadratic in the matter sources, S a b [see Eq. (0)]. Now, 

(5) 

we will add the contribution of E ab , which, as we have 
seen above, only contains the scalar part U . In order to 
study the dynamics of these models we will closely fol- 
low the analysis carried out in Paper I. This procedure 
is based on the introduction of new variables which de- 
scribe the dynamics in such a way that the state space is 
a compact space (see [ pT]Jl^| for details on the techniques 
of dynamical systems analysis). To carry out this study 
we have to consider four differentiated cases according to 
the signs of 3 R and U: (A) U > and 3 R < 0; (B) U > 
and 3 R > 0; (C) U < and 3 R < 0; (D) U < and 
3 R > 0. 



A. Case U > and i R < 

The state space in this case can be described by the 
dimensionless general-relativistic variables 



p ~ 3H 2 ' 



n k = 



>R 



6H 2 ' 

and the following new variables 
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Here fi p is the ordinary density parameter and 
Q\, and are the fractional contributions of the cur- 
vature, cosmological constant, brane tension, and the 
dark radiation energy, respectively, to the universe ex- 
pansion (^|). In fact, the Friedmann equation (|^) becomes 



n k + n A + n> 



Since these terms are non-negative they must belong 
to the interval [0, 1] and hence, the variables £1 = 
(f2p, f2fc, l^Aj Qu) define a compact state space. In- 
troducing the time derivative \H\~ l d/ dt, which we will 
denote by a prime, the evolution equation for H [Eq. (pi])] 
will decouple from the equations for SI. Then, the set of 
equations describing the dynamics is given by 



^ = e[2(l + g)-3 7 ]O p , 

SY k = 2eqn k , 

Q' A = 2e(l + q)Q A , 

Si' x = 2e(l + q - 3 7 )r> A , 

Q' u = 2e(q - l)Sl u , 



(15) 
(16) 
(17) 
(18) 
(19) 



III. DYNAMICS OF THE FLRW MODELS 



where q is 
q = 



3 7 -2 



n p - n A + (3 7 - i)n x + Mi 



In Paper I we studied the dynamics of the FLRW cos- 
mological models in a brane-world scenario in which the 



and e is the sign of H . As is clear, for e = 1 the model will 
be in expansion, and for e = — 1 it will be in contraction. 



3 



The critical points of the dynamical system (|15|)-(19), 
their coordinates in the state space and their eigenvalues 
are given in the following table |ll| 



Model 


Coordinates 


Eigenvalues 


F £ 


(1,0,0,0,0) 


e(3 7 - 2, 3 7 - 2, 37, -3 7 , 3 7 - 4) 


M c 


(0,1,0,0,0) 


e(-(3 7 - 2), 0,2,-2(3 7 -l), -2) 


dS E 


(0,0,1,0,0) 


-c(37,2,2,67,4) 


m £ 


(0,0,0,1,0) 


2e(^,(3 7 -l), 37,37-1,37 -2) 


R, 


(0,0,0,0,1) 


e(-(3 7 -4),2,4,-2(37-2),2) 



We have five hyperbolic critical points: the flat FLRW 
models (F), k = A = A" 1 = U = and a(t) = t 2 /^; 
the Milne universe (M), p = A= U = 0,k = — 1 and 
a(t) — t; the de Sitter model (dS), k = p = U = and 
a(t) — exp(-y/A/3 1); the non-general-relativistic BDL 
model (m) Ppp^| , found through a limiting process to be 
a(t) = i 1 ^ 3 ^ (see Paper I for more details); and a model 
(R) with the metric of a flat radiation FLRW model, 
a(t) = t 1 / 2 , but characterized by p = k = A = A -1 = 0. 



B. Case U > and i R > 



R e 



(e, 0,0, 1,0) 
(e, 0,0, 0,1) 



2e(37-l,^,37,37,3 7 -2) 
e(2, -(3 7 -4), 4, -2(3 7 -2), 4) 



Here E denotes a infinite set of saddle points whose line 
element is that of the Einstein universe (k = 1 and H — 
0). Their coordinates fl* satisfy (p0j) and 
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2fl* x ) - 2Q* U = 0. 



(21) 



In the table above 4> is given by 



37 
2 



(37-2)0; + 4 (37-1)0; 



4a* 



From the relation (f2l| ) and the allowed values of the vari- 
ables fi one can show that cj> is always positive for any 
7 . In GR the Einstein universe appears for 7 > 2/3 . In 
Paper I we showed that in brane-world scenarios without 
the contribution coming from the five-dimensional Weyl 
tensor it appears for 7 > 1/3 . Now, considering also this 
contribution Eq. ( f2l|) shows that the Einstein universe 
appears for any value of 7 . 



Here we can follow the study corresponding to the non- 
negative scalar curvature case in Paper I. Then, we will 
consider the variables tt = (Q,(l p , Cl\, ^w)> where Q 
is 



^ D ' 



D = H 2 



13 p 
6 K - 



and the variables with tilde are the analogues of those in 
( |i~3| ) and (O) but normalized with respect to D instead 
of H . Then, the Friedmann equation (Q) is now 



a, + n A + o A + n u = 1 



(20) 



Taking into account that q is now given by 

1 + qQ 2 = ^(O p + 20 A )+20 w , 

and using the time derivative ' = D~ 1 d/dt , the dynami- 
cal equations for f2 are decoupled from D and are given 
by 



Q' = -qQ 2 (l-Q 2 ), 
Cl' p = [2(1 + qQ 2 ) - 3 7 }QQ 
Q' A = 2(1 + qQ 2 )QCl A , 

n' x = 2[i + q Q 2 - 3-f]Qn x , 

h' u = -2(1 - qQ 2 )QVL u . 
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The critical points of this dynamical system, as well as 
their coordinates and eigenvalues are 



Model Coordinates 



Eigenvalues 



F e (e, 1,0, 0,0) e(3 7 - 2, 3 7 , 3 7 , -3 7 , 3 7 - 4) 

dS e (e, 0,1, 0,0) -e(2,3 7 ,0,67,4) 

e (o,n;,ni,ns;,n£) (o.v&o.-V&o) 



C. Case U < and 3 R < 

In this case we can get a compact state space by in- 
troducing the following set of dynamical variables Cl = 
(Z, Hp, Ofe, n A , where 



N 



N 



H 2 



2U 
A^ 



(22) 



and the variables with a bar are like those in ( |13| ) and ( |14| ) 
but normalized with respect to N instead of H . As be- 
fore, the Friedmann equation (JoJ) reduces to an expression 
in which all the terms are non-negative 



ftp + n k + n A + fix = 1 , 



(23) 



whence it follows that all of them are smaller than 1 . 
Using the time derivative ' = N~ 1 d/dt, in order to de- 
couple the evolution of N from the rest of variables, and 
taking into account the expression for q 

(g-l)Z 2 = f (ftp + 2ft A ) + ft fe -2, 

the dynamical system for f2 is 

Z' = -{q-l)Z 2 (\-Z 2 ) 1 

n' p = -i(3 1 -4)~2( q -i)z 2 }zn p , 

ft' fe = 2[1 + {q- l)Z 2 ]Zn k , 
Cl' A = 2[2 + (q - 1)Z 2 }ZQ A: 
Q'x = -2[(3 7 - 2) - (g - l)Z 2 ]ZClx • 

The critical points, their coordinates and eigenvalues are 
given in the table bellow fl3) 
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Model 


Coordinates 


F e 


(e, 1,0, 0,0) 


M £ 


(e, 0,1, 0,0) 


dS e 


(e, 0,0, 1,0) 


m c 


(e, 0,0, 0,1) 


S 


(0, Q* p , fl%, fi A , ^a) 



Eigenvalues 



e(3 7 - 4, 37,37-2,37,-37) 
e(-2, -(3 7 -2), 2, 2,-2(37- 1)) 
-e(4, 37,2,0,67) 
2e(3 7 -2,^,3 7 -l,37,3 7 ) 
(0, 0?, 0, 0, -v^?) 

We find new critical points, denoted by S, representing 
static models with vanishing and negative spatial curva- 
ture, i.e. H = and k = 0, -1 . Moreover fl* , fl* k , fl* A , 
and fl^ are real numbers satisfying ( p3|) and the following 
relation 



2fl* x ) = 0. 



The eigenvalues of S are given in terms of tp: 

^-¥[( 3 7-4)n; + 4(3 7 -2)nr 



20i. 



(24) 



(25) 



Eq. ( p4| ) implies that these static models (S) only appear 
for 7 > 2/3. Furthermore a, p and W are constants, 
(a* , p* ,U*), which must satisfy the following condition 



-U* > 



Xn 2 
2 



k 



(26) 



in order to have a positive energy density. 



D. Case U < and a R > 

This case is more involve since the Friedmann equa- 
tion (^J) has now two non-positive terms. Taking this 
into account we will consider the following dimensionless 
dynamical variables fl = (W, fl p , Oa, ^a> flu), where 



W 



H 
P 



P 2 = H 2 



2U 



fl p , fl\ , flu , and fl\ are the analogous of those in (|Tj|) 
and @ but normalized with P. It is important to note 
that flu does not appear in the Friedmann equation (|^), 
which now reads 



fi n + fi* 



fi, 



1 . 



(27) 



but it will be needed to close the dynamical system. 
Moreover, in this case flu is negative and from its defi- 
nition it belongs to the interval [—1,0] . Then, using the 
time derivative ' = P~ 1 d/dt, and the expression for q 

i + q w 2 = %(fi p + 2fi x ) + fi u , 

the evolution equations for fl are 

W = -[q(l - W 2 ) + fl u ]W 2 , (28) 
fl' p = [2(1 + qW 2 ) - (2fl u + &y)]w£l p , 

o A = 2(1 + q w 2 - n u )wn A , 

2[l + q w 2 - (O w + 3 7 )]^n A , 



fl' u = 2[qW 2 -flu- l]Wh u ■ 



The critical points and their coordinates and eigenvalues 
are 1131 



Model 


Coordinates 


Eigenvalues 




F e 


(e, 1,0, 0,0) 


e(3 7 - 2, 37, 37, -37, 37 


-4) 


dS e 


(e, 0,1, 0,0) 


-e(2,3 7 ,0,6 7 ,4) 




m £ 


(e, 0,0, 1,0) 


2 e (3 7 -2,^,37,3 7 ,37 


-2) 


E 


(o,n;,ni,nj,n^) 


(0,^,0,-^,0) 





where fl* , fl* A , fl x , and fl u are real numbers satisfying 
and 



2fl\ 



fit 



1 = 0. 



Now, the eigenvalues of the Einstein universe are deter- 
mined by ip, which is given by 



31 



i)fil 



4a* 



(3 7 -2)n; + 4(37- 

In contrast to what happens in the state space sector an- 



alyzed in subsection IIIB the Einstein universe now ap- 



pears for 7 > 1/3, just when ip starts being non-negative. 



E. Qualitative analysis 

Using the information obtained in the previous sub- 
sections we are going to describe qualitatively the full 
dynamics of FLRW in brane- world scenarios. To that 
end, we will construct and analyze the structure of the 
complete state space for these models. To start with, 
we need to know the dynamical character of the critical 
points. This can be deduced from the eigenvalues given 
above and summarized in the following table 



Model 




Dynamical character 






0< 7 < I 


7=1 


7>| 


F± 


saddle 


saddle 


saddle 


M± 


saddle 


saddle 


saddle 


dS+ 


attractor 


attractor 


attractor 


dS- 


repeller 


repeller 


repeller 


E 


saddle 


saddle 


saddle 


m+ 


saddle 


repeller 


repeller 


m_ 


saddle 


attractor 


attractor 


R+ 


repeller 


repeller 


saddle 


R- 


attractor 


attractor 


saddle 


S 




saddle 


saddle 



We can observe some changes with respect to GR and 
the analysis done for U = (Paper I). For instance, the 
Milne universe is always a saddle point for U ^ whereas 
M + (M_) is a repeller (attractor) for 7 < 2/3 in GR and 
for 7 < 1/3 in brane- world scenarios with U = . 

The full state space for the FLRW models is four- 
dimensional, due to the constraint imposed by the Fried- 
mann equa tion, and consists of four sectors (subsections 
III A - [II D| ) . They are disconnected in the sense that tra- 
jectories cannot change sector. Here, we will use a two- 
dimensional representation of the state space, formed by 
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G 




FIG. 3. S ector corr espon ding to the dynamics described in 



subsections [II A and [II B of the state space for the FLRW 
models with 7 £ (0, |) . It shows that the Einstein universe 
appears even for 7 < | , in contrast to what happens in GR 
and also when U = . 



all the two-dimensional invariant submanifolds, which in 
turn determine the whole dynamics. We have included 
the diagrams corresponding to the cases 7 6 (|, |) (see 
Fig. [I]) and 7 = | (see Fig. ||), which contain most of 
the physically interested cases, like dust (7=1) and ra- 
diation (7 = |). The second figure constitutes a new 
bifurcation^ due to the structure of the W-term, which 
is the same as that of a radiation fluid. Then, this bi- 
furcation is characterized by the presence of lines of new 
critical points with scale factor a(t) = t 1 / 2 but whose dy- 
namical character can be either that of a saddle, repeller 
or attractor, depending on the signs of 3 R and U. 

Coming back to the description o f the state space di- 

are 



agrams, the sectors of subsections [II B and III C 



represented by three rectangular areas, crossed by the 
dashed lines Q = and Z = respectively, enclosing 
tria ngular areas corresponding to the sector of subsec- 
tion [II A. The central area, crossed by the line W — 0, 



describes the dynamics of the sector in subsection [II D 



The shadowed region does not belong to the state space, 
it appears as a result of the two-dimensional representa- 
tion we have used. The dashed lines Q = Z = W = 0, 
where the Hubble function vanishes, separate regions of 



*The other bifurcations occur at 7 = 0, | (also for GR and 
U = 0) and 7 = i (also for 11 = 0). It is important to note 
that in the case 7 = | we find new lines of critical points 
corresponding to models with scale factor a(t) = t 1 ^ 2 . 



expanding models from regions with contracting models. 
Using this structure for the state space, the trajectories 
follow directly from the information given in the previous 
subsections. 

The attractors of the evolution are dS+ , R_ and m_ . 
More specifically, dS+ is an attractor for any value of 7, 
whereas R_ is an attractor only for 7 < 2/3 and m_ 
for 7 > 2/3. In fact, R_ and m_ exchange their role 
as an attractor at the bifurcation 7 = 2/3. Moreover, 
while dS + is an attractor for models that at some time 
will become ever expanding, R_ and m^ are attractors 
for (re)collapsing models. An interesting feature of the 
state space is that the volume of the region with mod- 
els evolving towards dS+ decreases as 7 increases, which 
means that the region of models collapsing in the future 
increases. The fact of having models that recollapse is 
directly connected with the appearance of static mod- 
els, namely E and S, which are saddle points. In this 
sense, it is important to remark that the presence of hi 
makes the Einstein models to appear for any value of 7 
(see Fig. ||), and the fact that it can be negative leads to 
the appearance of static models with non-positive spa- 
tial curvature, i.e. S. Furthermore, E and S are not 
isolated critical points but bidimensional surfaces in the 
state space. The location of these surfaces depending on 
the value of 7 has been represented in Fig. || 

The way of com pactif ying the s tate space for the sec- 
tors of subsections III A and III B is the same as the one 
used in Paper I for the sectors A R < and 3 R > re- 
spectively. The compactification of the sector in subsec- 

the 



tion 



III C is very similar to that of subsection III B 



only difference is that the negative term in the Friedmann 
equation is due to U instead of 3 R. In fact, as we can see 
in Figs, [l] and p, the structure is almost the same. For the 



sector corresponding to subsection [II D the situation is 



different due to the appearance of two non-positive con- 
tributions to the Friedmann equation. The compactifica- 
tion of the state space in this case leads to a new type of 
diagrams (see Fig. ||). To understand the dynamics we 
have to take into account that now we have the following 
invariant submanifold 



W 2 - ft 



u 



1 



(29) 



It divides the allowed rectangle [—1,1] x [—1,0] in the 
plane (W, ftu) into two regions with opposite sign of the 
spatial scalar curvature. From the definition of P we find 



W 2 - ftu - 1 = 



3 R 

Q P 2 



Since 3 R > 0, this means that only the interior region 
defined by ( |29| ) is allowed (see Fig. |5|). The special struc- 
ture of this sector of the state space leads to a new ex- 
ceptional feature consisting in the appearance of closed 
trajectories (see Figs, [l] and ||). They correspond to cos- 
mological models that after a certain period of time come 
back to the same situation as in the beginning of that 
period, and for this reason we will call them oscillating 



7 



n A A 



n A A 




FIG. 4. Surfaces of saddle critical points corresponding to static FLRW models (H = 0). The diagram on the left is the 
submanifold Z = of the sector 3 R < and U < and the critical points represented correspond to S. The diagrams on the 
center (right) show the location of the points E in the submanifolds W — (Q — 0) of the sector 3 R > and U < ( 3 R > and 
U > 0). Note that Q — 0, Z — and W — are not invariant submanifolds and therefore they do not contain any trajectories. 




m_ / m + 
tin ? fix yflu 

FIG. 5. Three-dimensional diagrams for the sector 3 _R > and U < of the state space of the FLRW models for 7 6 (|, |). 
Starting from the left, they correspond to the invariant submanifolds Cl\ = 0, Clu = and Cl\ — 0. The curved part of the 
boundary of the first and third state space diagrams correspond to the invariant submanifold W 2 — &u = 1. In these diagrams 
we also have closed trajectories describing oscillating universes. 
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universes. As we can see in Fig. [j] the trajectories fol- 
lowed by these models have a minimum and a maximum 
value of W which, from Eq. (plf), satisfy the following 
relation 



2fl x ) - U = -Si 



This equation has solutions only for 7 > i, just for the 
values of 7 that oscillating universes appear in the state 
space. Moreover, there is also a minimum and a max- 
imum of £).u that corresponds to W — 0. In fact, all 
the quantities describing these models are periodic and 
therefore have a minimum and a maximum value, in par- 
ticular the scale factor a(t) and the energy density p. 
Hence, these models do not have any spacelike singu- 
larity. In Fig. |^ we have represented the evolution of a 
particular class of oscillating universes. Apart from the 
features just described it is worth noting that while the 
transition from the expanding to the contracting stage is 
smooth, the converse is abrupt and mimics a non-singular 
bang or bounce. The avoidance of the singularity is due 
to the negative contribution of the dark-energy density 
term U, which goes like a -4 , to the universe expansion 
in the Friedmann equation (^). In the contracting stage 
this term grows and counterbalances the positive contri- 
butions in the Friedmann equation, stopping the expan- 
sion. This can only happen for ^ < 7 < | when the effect 
of the p 2 term, which goes like a -67 , is dominant, and 
for I < 7 < I (see Fig. when the effect of the p term, 
which goes like a~ 37 , dominates. After the bounce, the 
universe expands until the spatial curvature 3 i?, which 
goes like a -2 , stops the expansion to enter again a con- 
traction era. 



IV. DYNAMICS OF THE BIANCHI TYPE I 
MODELS 

In what follows we will complete the dynamical study 
of Bianchi type I models made in Paper I. Other works 
on these models in brane-world scenarios can be found 
L5l— GTt|| . Exact solutions on the brane have been con- 



sidered in |L8|], and a five-dimensional solution for the 
vacuum case has been presented in fl9| . 

As we have discussed above, the non-zero contributions 
from the five-dimensional Weyl tensor are 14 and V a b but 
since the second one has no evolution equation we will 
take V a b — 0. Then, we have to consider two differen- 
tiated cases according to the sign of U: (A) U > 0; (B) 
U < 0. 




p(t) 100 


40 

a 20- 









: ( 1 




I 1 



t (time) 

FIG. 6. Evolution of an oscillating FLRW universe with 
k = 1, 7 = 1 (dust), A = 0, and negligible contribution from 
the term quadratic in the energy density. We have taken k = 1 
and initial data: a(0) = H(0) = 1, p(0) = 9 and W(0) = —A. 



where we have defined a normalized variable for the shear 
contribution 



3H 2 



(30) 



The Friedmann equation @ then becomes 



tit 



1 . 



and again the state space is compact. Using the time 
derivative ' = \H\~ 1 d/dt the system of dynamical equa- 
tions is 



n; = e[2(l + g)-3 7 ]n />> 
Cl' A = 2e(l + q)Q A , 
n' a = 2e(q - 2)(l a , 
fi' A = 2e(l + q - 3 7 )r> A , 
fl' u = 2e(q - 1)Q U , 

where q has the following expression 
37-2. 



(31) 
(32) 
(33) 
(34) 
(35) 



q 



-n p + (3 7 - i)o A + 2Q CT - n A + n t 



A. Case U > 



The treatment of this case is very similar to that made 
in subsection III A. We will use the dimensionless vari- 
ables n = {n p ,n A ,n a ,n x ,n u ) [See E qs .© and @], 



The critical points of the dynamical system (|3l|)-(p5[), 
their coordinates in the state space and their eigenvalues 
are given in the next table 



Model 



Coordinates 



Eigenvalues 



F e 



(1, 0, 0, 0, 0) e(3 7 - 2, 3 7 , 3(7 - 2), -3 7 , 3 7 - 4) 
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dS e (0,1,0,0,0) 

K E (0,0,1,0,0) 

m £ (0,0,0,1,0) 

R e (0,0,0,0,1) 



-e(3 7 ,2,6,6 7 ,4) 
e(-3( 7 -2),6,4,-6(7-l),2) 
2e(^, 37, 3(7-1), 37 -1,37 -2) 
e(-(37-4),4,-2,-2(3 7 -2),2) 



The critical point K corresponds to the vacuum Kasncr 
models of GR, whose line element is given by (0) with 



A a =t 2p " and J2 Pa = J2 p * = 1 



(36) 



With respect to the analysis of Bianchi type I models 
done in Paper I, the only new critical points are the ra- 
diation models R± . 



B. Case U < 



The procedure we will follow in this case is the same 



as in subsection [II C . We will consider the following 

13) 



set of variables fi = (Z, Q pi Q\, Q ai Q\) [See Eq. 
where fl a is the analogous of Q a (^) but normalized with 
respect to N (p2|). With this variables the Friedmann 
equation (|j) reads 

n p + n A + n a + n x = 1 , (37) 

and the deceleration parameter q is given by 

(q - \)Z 2 = ?2(n p + 2fi A ) + 3ti a - 2 . 

With respect to the time derivative ' = N~ 1 d/dt, the 
dynamical system for ti is 



Z' 
9' 



-(q-l)Z 2 (l-Z 2 ), 
.„ ~[(3j-4)-2( q -l)Z 2 ]Zn p , 

n' A = 2[2 + { q -i)z 2 ]zn A , 
n' a = -2[i - (q - i)z 2 ]zn a , 

Cl' x = -2[(3 1 -2)-{q-l)Z 2 ]ZCl x . 



(38) 
(39) 
(40) 
(41) 
(42) 



The critical points of the dynamical system (B^)-(42), 
their coordinates in the stat e sp ace and their eigenvalues 
are given in the next table p3[ 



Model 



Coordinates 



Eigenvalues 



F e (e, 1,0, 0,0) 

dS e (e, 0,1, 0,0) 

K £ (e,0,0,l,0) 
m E (e,0,0,0,l) 

nE (0,fi;,Ql,Q%,Ql) 



e(3 7 -4,37,3 7 ,3( 7 - 2),-3 7 ) 

-6(4,37,0,6,67) 
£ (2, -3(7 -2), 6, 6, -6(7-1)) 
2e(37-2,^,3 7 ,3(7-l),3 7 ) 



(0, VC,o,o,-v/c) 

We have a new set of saddle critical points, designate d by 
nE, whose coordinates in the state space n* satisfy (|37j) 
and 

%(n* p + 2ni) + 3n* a -2 = o. 

The quantity £ is given by 



•4)n; + 4(3 7 - 



2)nx 



and in the FLRW limit, a = , it corresponds to the 
quantity ip in ( p5| ) for k — 0. The nE models are non- 
expanding (H = 0) Bianchi models which in the FLRW 
limit become static, in fact, in this limit, they coincide 
with the critical points S with 3 R — 0. A noteworthy 
difference is that the models S only appear for 7 > 2/3 
whereas the models nE appear for all values of 7. More- 
over, p, a and U are constants, (p*,cr*,U*), which, in 
order to have a positive energy density must satisfy a 
condition analogous to ( p6| ) 

6 V 

In fact, integrating the effective Einstein equations (|^), 
we can find the metric functions A a (t) in (^) correspond- 
ing to the critical points nE, 

A a {t) = e«- t , 

where q a (a = 1 — 3) are constants (proportional to the 
components of the shear tensor) such that 



3 

qa — and 



3 

a = l 



2a* 



Note that despite the scalar factor remains constant 
(]l a =i A <x(t) = !) the matter can expand or contract 
along the shear principal axes. In fact, when only one q a 
is negative the system evolves towards a pancake singu- 
larity, whereas when two of them are negative it evolves 
towards a cigar singularity (for details on this type of 
singularities see, e.g., CO]). 



C. Qualitative analysis 

To analyze the dynamics we have to consider first the 
character of the critical points, which can be deduced 
from the information given in the previous subsections 
and can be summarized in the following table 



Model 




Dynamical character 






< 7 < 1 


7 = 1 


7 > 1 


F± 


saddle 


saddle 


saddle 


dS+ 


attractor 


attractor 


attractor 


dS_ 


repeller 


repeller 


repeller 


K+ 


repeller 


repeller 


saddle 


K_ 


attractor 


attractor 


saddle 


m+ 


saddle 


repeller 


repeller 


m_ 


saddle 


attractor 


attractor 


R± 


saddle 


saddle 


saddle 


nE 


saddle 


saddle 


saddle 



The character of the critical points F±, dS±, K±, and 
m± is the same as in the case U = studied in Paper I. 
The new dynamical features are due to the appearance 
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of the critical points R± and nE. Moreover, we also find 
two new bifurcations at 7 = | and 7 = | , besides the 
bifurcations at 7 = 0,1,2 already found in Paper I for 
the case hi = 0. 

The state space for Bianchi type I models is also four- 
dimensional and is divided into two sectors according to 
the sign of hi. We have used the same type of bidimcn- 
sional representation in which only the trajectories for 
the two-dimensional invariant submanifolds are drawn. 
In Fig. ^ we show the state space for 7 6 (1, |), which 
includes physically well-motivated situations from dust to 
radiation universes. The sector hi < described in sub- 



section [VB is represented by three rectangular regions 



which enclose a tria ngle w here the dynamics of the sector 
hi > (subsection IV A ) is depicted. The dashed lines 
Z = divided regions of expanding models from regions 
of contracting models, and these are the places where 
the non-expanding models nE appear. In fact, these sad- 
dle points are not isolated points but two-dimensional 
surfaces. The location of the surfaces of points nE for 
different values of 7 has been drawn in Fig. 0. 

For the sector hi > the general attractor is de Sit- 
ter, as it happens when hi = 0. The situation changes 
for the sector hi < 0, where the possible final states of 
the evolution are dS+ , K_ and m_ . Like in some sectors 
of state space of the FLRW models, the sector hi < 
of Bianchi models is divided into two regions, the re- 
gion of models that will approach dS+ and then will 
become ever expanding, and the region of models that 
will (re)collapse. Obviously the models approaching dS+ 
in the future isotropize, although there is an interme- 
diate stage in which the relative contribution of shear 
grows and reaches a maximum. In the sector hi > this 
maximum is reached when the following condition holds 
(see0) 



ficr — 1 — ^flll 



and in the sector hi < the condition is 



a, 



1 



2n> 



On the other hand, the attractor for collapsing models 
(hi < 0) is in general K (m_ for the submanifold a = 0) 
for < 7 < 1 and m_ (K_ for the submanifolds ti\ = 0) 
for 7 > 1. For dust (7 = 1) we have a bifurcacion in 
which the attractors are anisotropic models whose line 
element is the same as the one of Kasner models [see 
Eqs.(f7|),(|36l)l but the constant p a now only satisfy 



1 . 



From this discussion we deduce that there are sets of tra- 
jectories in the state space that either do not isotropize 
in the future or isotropize but the final state is not de 
Sitter (dS+). Hence, the cosmic no- hair theorem is not 
satisfied in the hi < sector. This fact has been an- 
ticipated in the analysis made in |jq| , although in their 




FIG. 7. Complete state space for Bianchi type I models 
with 7 G fl. 4 ' 



4> i 

'-> 3' 



non-compact representation of the state space the final 
state for collapsing models was not shown. 

Finally, with regard to the structure of the initial sin- 
gularity, in Paper I we showed that for hi — and 7 > 1 
the initial singularity is isotropic, in contrast to what 
happens in GR. Now, this is also true for the sector 
hi > 0. In the sector hi < the singularity is also isotropic 
except for the submanifold Cl\ = 0. 



V. REMARKS AND CONCLUSIONS 

In this paper we have completed the dynamical sys- 
tems analysis on the evolution of cosmological models in 
brane-world type scenarios started in Paper I, where the 
effects of the five-dimensional (bulk) Weyl tensor were 
not considered. By including this quantity in the anal- 
ysis we have constructed the most complete state space 
for non-tilted perfect-fluid FLRW cosmological models. 
In the case of Bianchi type I models the state space we 
have presented is also complete except for the fact that we 
have disregarded the components of the five-dimensional 
Weyl tensor for which the theory does not provide evolu- 
tion equations, namely V a b- This means that we cannot 
make predictions on how the term V a b affects the dy- 
namics. Perhaps future studies on string/M-theory could 
shed light on this question. 

From the dynamical study carried out we have found 
new interesting features not present in GR nor in brane- 
world scenarios with hi = 0. In particular, we have dis- 
covered new critical points: R±, S and nE. The Einstein 
universe E and the non-expanding anisotropic models nE 
appear for any value of 7, which is directly connected to 



11 




brane-world scenarios the four-dimensional cosmological 
constant is given by 



FIG. 8. Representation of surfaces of saddle critical points 
corresponding to the non-expanding models nE in the Z = 
submanifold of the state space of Bianchi type I models. 



the fact that we have now (re)collapsing models both for 
FLRW and Bianchi type I models for any 7. It is also 
worth noting that now we have static critical points for 
any sign of the spatial curvature 3 R. We have also seen 
that de Sitter (dS+) is always an attractor but not the 
only one. For FLRW models, R_ and m_ are also attrac- 
tors for 7 < I and 7 > | respectively, and for Bianchi 
type I models, K_ and m_ are attractors for 7 < 1 and 
7 > 1 respectively. 

A relevant interesting feature is the appearance of os- 
cillating universes in the 3 R > and U < sector of 
the state space of the FLRW models. In these universes 
the physical variables oscillate periodically between a 
minimum and a maximum value without reaching any 
spacelike singularity. With regard to the evolution of the 
anisotropy in Bianchi type I models we have found that 
for U > they always isotropize whereas for U < they 
can both isotropize or collapse. This violation of the cos- 
mic no-hair theorem has also been pointed out in [ jl6| . 
The use in the present paper of dynamical variables that 
compactify the state space allows us to see clearly which 
are the possible final states of the collapse (K_ or m_) 
and how the state space is divided into two regions ac- 
cording to whether the models isotropize or collapse. On 
the other hand, we have seen that the initial singularity 
is isotropic for 7 > 1, as it also happens in the particular 
case U = 0. 

Finally, the analysis we have carried out in this pa- 
per provides information on the question of the stability 
of the so-called Randall- Sundrum fine-tuning condition, 
which consists in the vanishing of the four-dimensional 
cosmological constant (A = 0). More specifically, in 



A 



|A (5 



Xi = 6 



|A(5)| 
K (5) 



where A c is a critical tension. Then the fine tuning con- 
sists in adjusting the brane tension A so that A = A c . 
From our perspective, we can deal with this question by 
asking whether or not the invariant submanifold JIa = 
(and also the equivalent ones fl\ = = = 0), where 
the fine-tuning condition holds, is dynamically stable. In 
other words, what happens when initial conditions are 
prescribed near this invariant submanifold of the state 
space. It can be seen that for the regions of recollapsing 
models the evolution will bring them back to f?A = 0, 
although if we start with H > the models will initially 
separate from that submanifold. For the other regions, 
the evolution will move the models away from fl\ = 
in such a way that they will approach de Sitter, which 
is precisely the opposite situation to A = (see also plf] 
for an alternative discussion). Then, we can conclude 
that the fine-tuning condition is not dynamically stable 
in the sense that there are non-zero measure regions of 
initial data near the fine-tuning A = evolving towards 
a situation in which that condition does not longer hold. 



Acknowledgements: The authors wish to thank their 
colleagues in the Relativity and Cosmology Group for 
fruitful discussions. This work has been supported by the 
European Commission (contracts HPMF-CT-1999-00149 
and HPMF-CT-1999-00158). 



[1] L. Randall and R. Sundrum, Phys. Rev. Lett 83, 4690 
(1999). 

[2] P. Hofava and E. Witten, Nucl. Phys. B460, 506 (1996); 
P. Hofava and E. Witten, Nucl. Phys. B475, 94 (1996). 

[3] K. Benakli, Int. J. Mod. Phys. D 8, 153 (1999); A. Lukas, 
B. A. Ovrut, K. S. Stelle, and D. Waldram, Phys. Rev. D 
59, 086001 (1999); H. S. Reall, Phys. Rev. D 59, 103506 
(1999); A. Lukas, B. A. Ovrut, and D. Waldram, Phys. 
Rev. D 60, 086001 (1999); A. Lukas, B. A. Ovrut, and 
D. Waldram, Phys. Rev. D 61, 023506 (1999); A. P. Bill- 
yard, A. A. Coley, J. E. Lidsey, and U. S. Nilsson, Phys. 
Rev. D 61, 043504 (2000). 

[4] P. Binetruy, C. Deffayet, and D. Langlois, Nucl. Phys. 
B565, 269 (2000). 

[5] P. Binetruy, C. Deffayet, U. Ellwanger, and D. Langlois, 
Phys. Lett. B 477, 285 (2000). 

[6] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D 
62, 024012 (2000). 

[7] M. Sasaki, T. Shiromizu, and K. Maeda, Phys. Rev. D 
62, 024008 (2000). 



12 



[8] A. Campos and C. F. Sopuerta, Phys. Rev. D 63, 104012 
(2001). 

[9] R. Maartens, Phys. Rev. D 62, 084023 (2000). 
[10] S. W. Hawking and G. F. R. Ellis, The large scale struc- 
ture of space-time (Cambridge University Press, Cam- 
bridge, 1973). 

[11] J. Wainwright and G. F. R. Ellis, Dynamical systems 
in cosmology (Cambridge University Press, Cambridge, 
1997); A. A. Coley, in Proceedings of the Spanish Relativ- 
ity Meeting. ERE-99, edited by J. Ibafiez (Euskal Herriko 
Unibertsitatea, Bilbo, 2000). 

[12] D. K. Arrowsmith and C. M. Place, An introduction to 
dynamical systems (Cambridge University Press, Cam- 
bridge, 1990); D. W. Jordan and P. Smith, Nonlinear or- 
dinary differential equations. An introduction to dynam- 
ical systems (Oxford University Press, Oxford, 1999). 

[13] The presence of a zero eigenvalue is (excepting for bi- 
furcations and the sets of critical points E, S and nE) a 
consequence of having a dynamical system with a con- 
straint: The Friedmann equation. If we had used it to 
eliminate one variable and reduce the dimension of the 
state space, we would not find any vanishing eigenvalue. 
In conclusion, this fact does not affect to the hyperbolic 
character of the critical points. 

[14] E. E. Flanagan, S.-H. H. Tye, and I. Wasserman, Phys. 
Rev. D 62, 044039 (2000). 

[15] R. Maartens, V. Sahni and T. D. Saini, Phys. Rev. D 63, 
063509 (2001). 

[16] M. G. Sant os, F. Vernizzi, and P. G. Ferreira, jiep 



ph/0103112. 



[17] A. V. Toporensky, gr-qc/0103093 



[18] C.-M. Chen, T. Harko, and M . K. Mak, ^iep-th/010324[ . 
[19] A. V. Frolov, |gr-qc/0102064| 

[20] H. Stephani, General Relativity (Cambridge University 

Press, Cambridge, 1990). 
[21] T. Boehm, R. Durrer, and C. van de Bruck, |iep 



th/0102144 



13 



